Experimental reconstruction of photon statistics without photon counting 
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Experimental reconstructions of photon number distributions of both continuous- wave and pulsed 
light beams are reported. Our scheme is based on on/off avalanche photodetection assisted by 
maximum-likelihood estimation and does not involve photon counting. Reconstructions of the dis- 
tribution for both semiclassical and quantum states of light are reported for single-mode as well as 
for multimode beams. 



PACS numbers: 42.50.Ar, 42.50.Dv, 03.65.Wj 

The measurement of the statistical distribution of the 
number of photons provides fundamental information on 
the nature of any optical field. The choice of a detec- 
tor with internal gain suitable for the measurement is 
not trivial when the flux of the photons to be counted is 
such that more than one photon is detected in the time- 
window of the measurement, which is set by the detector 
pulse-response, or by an electronic gate on the detector 
output, or by the duration of the light pulse. In this case, 
we need a congruous linearity in the internal current am- 
plification process: each of the single electrons produced 
by the different photons in the primary step of the detec- 
tion process (either ionization or promotion to a conduc- 
tion band) must experience the same average gain and 
this gain must have sufficiently low spread. The fulfill- 
ment of both requisites is necessary for the charge integral 
of the output current pulse be proportional to the number 
of detected photons. Photon detectors that can operate 
as photon counters are rather rare. Among these, Photo- 
Multiplier Tubes (PMT's) and hybrid photodetectors 
have the drawback of a low quantum efficiency, since 
the detection starts with the emission of an electron from 
the photocathode. Solid state detectors with internal 
gain, in which the nature of the primary detection pro- 
cess ensures higher efficiency, are still under development. 
Highly efficient thermal detectors have also been used as 
photon counters, though their operating conditions are 
still extreme (cryogenic conditions) to allow common use 
0, Q . The advent of quantum tomography provided an 
alternative method to measure photon number distribu- 
tions 5j. However, the tomography of a state, which 
has been applied to several quantum states [|| , needs the 
implementation of homodyne detection, which in turn re- 



quires the appropriate mode matching of the signal with 
a suitable local oscillator at a beam splitter. Such mode 
matching is a particularly challenging task in the case of 
pulsed optical fields. 

Photodetectors that are usually employed in quantum 
optics such as Avalanche PhotoDiodes (APD's) operat- 
ing in the Geiger mode 0, Q seem to be by definition 
useless as photon counters. They are the solid state pho- 
todetectors with the highest quantum efficiency but they 
have the obvious drawback that the breakdown current 
is independent of the number of detected photons, which 
in turn cannot be determined. The outcome of these 
APD's is either "off" (no photons detected) or "on" i.e. 
a "click", indicating the detection of one or more pho- 
tons. Actually, such an outcome can be provided by 
any photodetector (PMT, hybrid photodetector, cryo- 
genic thermal detector) for which the charge contained 
in dark pulses is definitely below that of the output cur- 
rent pulses corresponding to the detection of at least one 
photon. Note that for most high-gain PMT's the anodic 
pulses corresponding to no photons detected can be easily 
discriminated by a threshold from those corresponding to 
the detection of one or more photons. 

The statistics of the " no-click" and " click" events from 
an on/off detector, assuming no dark counts, is given by 



P0(V) = 5>-77rVn 



(1) 
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and p>o{rf) = 1 — Poiv) i where g n is the probability of 
finding n photons and r\ is the quantum efficiency of the 
detector, i.e. the probability of a single photon to be 
revealed. At first sight the statistics of an on/off detec- 
tor appears to provide quite a scarce piece of informa- 
tion about the state under investigation. However, if the 
statistics about po(r)) is collected for a suitably large set 
of efficiency values then the information is enough to re- 
construct the whole photon distribution g n of the signal, 
upon a suitable truncation of the Hilbert space. 
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The reconstruction of photon distribution through 
on/off detection at different efficiencies has been analyzed 
@ and its statistical reliability investigated in some de- 
tails . In addition, the case of few and small values of r\ 
[To| has been addressed. However, whilst these theoreti- 
cal studies found an application to realize a multichannel 
fiber loop detector |lllll2| . an experimental implementa- 
tion of this technique for reconstructing photon distribu- 
tion of a free-propagating field is still missing. In view of 
the relevance of photon distribution for applications in 
quantum information and foundations of quantum me- 
chanics, the purpose of this letter is to show that a recon- 
struction of the photon distribution by using this tech- 
nique can be effectively realized experimentally with ex- 
cellent results both for free-propagating continuous- wave 
(cw) and pulsed light beams, for both single-mode semi- 
classical and quantum states, as well as for multimode 
states. 

The procedure consists in measuring a given signal by 
on/off detection using different values r) v (v = 1, ...,K) 
of the quantum efficiency. The information provided by 
experimental data is contained in the collection of fre- 
quencies f v = fo(r) v ) = TiQv/nv where no is the number 
of "no click" events and n v the total number of runs with 
quantum efficiency r\ v . Then we consider expression £Q) 
as a statistical model for the parameters g n to be solved 
by maximum-likelihood (ML) estimation. Upon defining 
p v = pa{rj v ) and A vn — (1 — r} v ) n we rewrite expres- 
sion Q as p v — '^2 in A vn g n . Since the model is linear 
and the parameters to be estimates are positive (LINPOS 
problem) , then the solution can be obtained by using the 
Expectation-Maximization algorithm (EM) [jjj. By im- 
posing the restriction J2 n Qn = 1> we obtain the iterative 
solution 
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where Pi/[{gn^}] are the probabilities p v , as calculated 
by using the reconstructed distribution {g^ } at the 
j-th iteration. As a measure of convergence we use 
the total absolute error at the i-th iteration gW = 
S^=o ~ P»[{#n }] and stop the algorithm as soon as 
eW goes below a given level. The total error measures 
the distance of the probabilities Pi/[{f?n^}], as calculated 
at the i-th iteration, from the actual experimental fre- 
quencies. As a measure of accuracy we adopt the fidelity 

q{i) — J2 n \J g n q$ between the reconstructed distribu- 
tion and the theoretical one. 

In order to verify the potentialities of this technique 
we applied it to the reconstruction of various quantum 
optical states, generated either in the cw or in the pulsed 
regimes. As a first example we have considered single 
photon states that have been generated by producing 
Parametric Down Conversion (PDC) heralded photons. 
In little more detail, a pair of correlated photons has 
been generated by pumping a type II /3-Barium-Borate 
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FIG. 1: Reconstruction of the photon distribution for a weak 
coherent state (grey) and for the heralded single-photon state 
produced in type II PDC (black). Inset: Experimental fre- 
quencies /„ of no-click events as a function of the quantum 
efficiency r\ v for a weak coherent state (upper curve) and for 
PDC heralded photon state (lower curve) compared with the 



theoretical curves, p„ 
tively. 



I — r\ u \a}\ and p v = 1 — r\ v respec- 



(BBO) crystal with a cw Argon ion laser beam (351 nm) 
in collinear geometry. After having split the photons of 
the pair by means of a polarizing beam splitter, the detec- 
tion of one of the two by a silicon avalanche photodiodc 
detector (SPCM-AQR-f 5, Perkin Elmer) was used as an 
indication of the presence of the second photon in the 
other channel, namely a window of 4.9 ns was opened 
for detection in arm 2 in correspondence to the detec- 
tion of a photon in arm I. This "heralded photon" was 
then measured by a silicon avalanche photodiode detec- 
tor (SPCM-AQR-I5, Perkin Elmer) preceded by an iris 
and an interference filter (IF) at 702 nm, 4 nm FWHM, 
inserted with the purpose of reducing the stray light. The 
quantum efficiency of the detection apparatus (including 
IF and iris) was measured to be 20% by using the PDC 
calibration scheme (see |14|'). Lower quantum efficiencies 
were simulated by inserting calibrated neutral optical fil- 
ters on the optical path. A comparison of the observed 
frequencies /„ with the theoretical curve (1 — rj u ) is pre- 
sented in the inset of Fig. The photon distribution 
has been reconstructed using K = 34 different values of 
the quantum efficiency from r\ v ~ to r\ v ~ 20% with 
n v = I0 6 runs for each r\ u . Results at iteration i = 10 6 
are shown in Fig. ^ As expected, the PDC heralded 
photon state largely agrees with a single photon Fock 
state. However, also a small two photons component 
and a vacuum one are observed. The p2 contribution 
is expected, by estimating the probability that a sec- 
ond photon randomly enters the detection window, to 
be 1.85% of pi, in agreement with what observed. A non 
zero po is also expected due to background. This quan- 
tity can be evaluated to correspond to (2.7 ± 0.2)% by 
measuring the counts when the polarization of the pump 
beam is rotated to avoid generation of parametric fluo- 
rescence. Also this estimate is in good agreement with 
the reconstructed po. A second example is represented 
by a strongly attenuated coherent state, which has been 
produced by a He-Ne laser beam attenuated to photon- 
counting regime by insertion of neutral filters. Also in 
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this case the same silicon avalanche photodiode detector 
was used. The reconstructed distribution (with K = 15 
different values of the quantum efficiency from 7/„ ~ to 
r\ v ~ 66% with n v — 10 6 runs for each rj v ) agrees well 
with what expected for a coherent state with average 
number of photons \a\ 2 ~ 0.02. In the inset of Fig. Q]the 
frequencies f v as a function of r\ v are compared with the 
theoretical prediction p u = exp{— rj v \a\ 2 } ~ 1 — i] v \a\ 2 . 
Notice that in this case we do not have IF or irises in 
front of the detector and all the other attenuations can 
be included in the generation of the state: thus the high- 
est quantum efficiency is taken to be 66% as declared by 
the manufacturer data-sheet for the photodetector. 

In the pulsed domain, we have measured three different 
optical states generated starting from the third harmon- 
ics (349 em, 4.45 ps) of a cw mode-locked Nd:YLF laser 
regeneratively amplified at a repetition rate of 500 Hz 
(High Q Laser). For all the measurements, the light was 
delivered to a photo multiplier tube (PMT, Burle 8850) 
through a multi mode fiber (100 fim core diameter). Al- 
though the PMT has the capability of counting the num- 
ber of photoelectrons produced by one or more photons 
0, for the present application we used it in a Geiger- 
like configuration, by setting a threshold to discriminate 
on/off events. Furthermore, by using the detector in the 
regime of linear response, the knowledge of its photocath- 
ode quantum efficiency is sufficient ti determine the rjv 
values. 

The first measurement was performed on the pulse 
emerging from the laser source. Due to the pulsed nature 
of the source, we do not expect to recover a true Poisso- 
nian statistics. Rather, we expect a Gaussian distribu- 
tion g n ,G IQ5 with mean value N and variance N + a 2 
which takes into account the presence of noise. N is the 
photon mean value and a 2 /N is as a measure of the devi- 
ation from Poissonian statistics. In Fig. [3] a) we show the 
photon distribution g n , reconstructed at the i = 50000 
iteration of the ML algorithm, along with the best fit ob- 
tained with the model g n> o (fitting parameters N = 4.88 
and a 2 = 0.63). The inset of the figure compares the 
experimental frequency /„ data (K — 37 values of 77, 
n v = 10 4 runs for each rf) as a function of rj v with the 
theoretical values calculated through and the param- 
eters given by the fit of the photon distribution. Both 
the reconstructed distribution and the experimental fre- 
quencies agrees very well with the above Gaussian model. 
The fidelity of the reconstruction is G ~ 0.998. Using the 
estimated value of a 2 , a deviation of about 13% from the 
Poissonian statistics can be derived for the laser photon 
number distribution. 

A second measurement was performed on laser pulses 
diffused by a moving ground glass. If the photons are 
collected from within an area of spatial coherence, the 
system acts as a pseudo-thermal source, whose photon 
number distribution is given by g n> x = N n (N + l) _,l_1 . 
Figure [21 b) shows the photon distribution g n , as recon- 
structed at the 400th iteration and the best fit of the data 
with g n: T (N = 5.33); the fidelity is given by G ~ 0.995. 
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FIG. 2: Reconstructed photon distribution (black bars) and 
best theoretical fit (grey bars) for three different states in the 
pulsed regime: a) laser pulse, b) diffused laser pulse, c) mul- 
timode state produced in type I PDC. Insets: Experimental 
frequency /„ data in function of the quantum efficiency r\ v 
and theoretical model for each one of the states. 



The inset of the figure contains the experimental fre- 
quency f v data [K = 24 values of 77, n„ = 10 4 runs for 
each 77) and their theoretical values as calculated from 

CD- 

The last measurement was performed on the blue por- 
tion (420 nm) of the down conversion fluorescence pro- 
duced by a type I BBO crystal (10 mm depth, cut at 
34 deg) pumped by the laser pulse. The pump, incident 
orthogonally to the crystal face, had an intensity ~ 60 
GW/cm . In this experimental condition we expect a 
coherence time of the generated field of ~ 1 ps, that 
corresponds to measuring a convolution of 4-5 temporal 
modes I lot . The photon number distribution is expected 
to be a multithcrmal distribution of the form 

= (n + M-1)! ,„s 

9n ' M n\(ji- 1)1(1 + N/n)»(l+n/N)" ' ( > 

where /j is the number of temporal modes. The photon 
distribution reconstructed at the i = 1500 iteration, is 
shown in Fig. [5] c) along with the best fit of the data 
using © (N = 6.17 and fx = 5); the fidelity of the recon- 
struction is given by G ~ 1. In the inset of the figure we 
show the experimental frequency /„ data (K = 18 val- 
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ues of r\, n„ = 10 4 runs for each 77) and their theoretical 
values as calculated according to Q . 

As a comment to the experimental results in the pulsed 
regime, we note that the best reconstruction of the pho- 
ton distribution is achieved at a different number of it- 
erations for the three different measured optical states, 
and that the absolute error e does not approach the same 
value. This is due to the presence of excess noise in our 
measurements, since the stability and the repetition rate 
of our source (500 Hz) limits to n v ~ 10 4 the number 
of runs for each value of the quantum efficiency @ . The 
choice of the best iteration to stop the algorithm is driven 
by the possibility to fit the distribution with a suitable 
model. We stress that there was no a-priori decision in 
choosing a Gaussian distribution for case a) or of a multi- 
thermal distribution for case c), but, on the contrary, we 
followed the a-posteriori observation that no other dis- 
tribution could fit equally well the reconstructed data. 

As to the comparison of the present technique with 
other schemes to reconstruct the photon distribution we 
have to distinguish between cw and pulsed regime. For 
cw field an alternative technique is represented by quan- 
tum homodyne tomography (QHT), which has indeed 
been applied to the reconstruction of single-photon and 
single-photon added states of light 0] . The advantage of 
our technique compared to QHT is twofold. On one hand 
QHT requires, for the same task, a more complex appa- 
ratus and high-efficiency photodetectors. On the other 
hand, QHT is more noisy. In fact, homodyne data con- 
tains the whole information about the state under inves- 



tigation (not only the photon distribution) and, in turn, 
this results in a more noisy determination when only part 
of the information is of interest . In order to obtain 
results such those of Fig. Q]QHT requires a by far larger 
set of data. In the pulsed regime, where realization of 
QHT is still challenging, the direct measurement of the 
photon statistics can be conveniently done by PMT's or 
hybrid photodetectors [l|,l2|. However, owing to limita- 
tions of the values of the maximum detection efficiency, 
results at number of photons as low as those in Fig. 2 
cannot be obtained by PMT's and are at limit of fea- 
sibility for hybrid photodetectors. While our method is 
particularly suited to measure the photon distributions of 
low-intensity field, where the above detectors are not ef- 
fective, it can be always applied to higher intensity fields. 

In conclusion, we experimentally implemented a recon- 
struction method for the photon distribution based on 
on/off detection at different quantum efficiency followed 
by a ML iterative algorithm. Our experimental results 
demonstrate that the technique can be applied to both 
cw and pulsed regimes, and for a wide range of signal en- 
ergy (from single- photon states to mesoscopic signals), a 
feature that makes it preferable to other methods till now 
devised for reconstructing photon number distribution. 
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